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Abstract. Recently, several complex network approaches to time series analysis have been developed 
and applied to study a wide range of model systems as well as real- world data, e.g., geophysical or fi- 
nancial time series. Among these techniques, recurrence-based concepts and prominently e-recurrence 
networks, most faithfully represent the geometrical fine structure of the attractors underlying chaotic 
(and less interestingly non-chaotic) time series. In this paper we demonstrate that the well known graph 
theoretical properties local clustering coefficient and global (network) transitivity can meaningfully be 
exploited to define two new local and two new global measures of dimension in phase space: local upper 
and lower clustering dimension as well as global upper and lower transitivity dimension. Rigorous an- 
alytical as well as numerical results for self-similar sets and simple chaotic model systems suggest that 
these measures are well-behaved in most non-pathological situations and that they can be estimated 
reasonably well using e-recurrence networks constructed from relatively short time series. Moreover, we 
study the relationship between clustering and transitivity dimensions on the one hand, and traditional 
measures like pointwise dimension or local Lyapunov dimension on the other hand. We also provide 
further evidence that the local clustering coefficients, or equivalently the local clustering dimensions, 
are useful for identifying unstable periodic orbits and other dynamically invariant objects from time se- 
ries. Our results demonstrate that e-recurrence networks exhibit an important link between dynamical 
systems and graph theory. 



I Introduction 

Recurrence of previous states is a key property of dy- 
namical systems [T]. In mathematical terms, one speaks 
of a recurrence if at time tj, a trajectory of a complex 
system returns into the dynamical neighbourhood of 
a previous state Xi = x{U) (ti < tj). Such a neighbour- 
hood can be defined by either considering the /c-nearest 
neighbours of Xi (fixed mass of the neighbourhood) or 
in terms of an e-ball B^{xi) centered at Xi (fixed volume 
of the neighbourhood). In the latter case, one defines 
the binary recurrence matrix 

R,,{e)^0{e^\\x,-x,\\) (1) 

for a trajectory sampled at a fixed number of points 
ti in time, where Q{-) is the Heaviside function and 

II • II some norm (e.g., maximum or Euclidean norm) 
in the metric space including the considered attractor. 
From this definition, it is evident that a recurrence di- 
rectly corresponds to i?y — 1. For the sake of clarity, 
we will more specifically speak about an e-recurrence 
in the following to distinguish this definition from the 
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alternative fc-nearest neighbour based definition. The 
recurrence matrix can be easily visualised in terms of 
recurrence plots (RPs) II1I2I3L the appearance of which 
allows for a simple graphical discrimination between 
qualitatively different types of dynamics. 

Since the poineering work of Poincare [4| on the 
three-body problem, it has been more and more recog- 
nised that recurrences are important for understand- 
ing the overall dynamics [3|, i.e., they encode all rele- 
vant information about the behaviour of a stationary 
system. Even more, it has been proven recently that 
the temporal pattern of recurrences allows us to recon- 
struct the rank-order of any scalar time series (i.e., a 
unique representation of the original t rajecto ry up to a 
scaling by some monotonous function) 1516171 . As a con- 
sequence, many well established dynamical invariants 
(such as correlation entropy, correlation dimension, or 
2nd-order mutual information) can be estimated from 
RPs [8 |. Moreover, the statistical analysis of line struc- 
tures in RPs known as recurrence quantification anal- 
ysis (RQA) allows the definition of a large variety of 
additional measures of complexity, which correspond 
to time intervals with similar consecutive states (ver- 
tical/horizontal lines) or time evolution (diagonal lines 
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off the main diagonal), respectively. Since these mea- 
sures are easily calculable, RQA has found numerous 
applications in the last two decades [1|. 

Besides the phenomenon of recurrence, another ap- 
pealing paradigm that has attracted considerable in- 
terest over the last years are networks with complex 
topology, called shortly complex networks. In the most 
simple case (to which we will restrict our attention in 
this work), a network can be mathematically described 
as simple graph G = {V, E), where V = {!,..., N} is 
the set of vertices with \V\ = N, and EC V x V is 
the set of edges between pairs of vertices. Note that we 
do neither consider multiple edges between two ver- 
tices, nor hyperedges connecting more than two ver- 
tices with each other, i.e. \E\ < N{N - l)/2. Further- 
more, we will restrict our attention to unweighted and 
undirected networks. In this case, the whole network 
connectivity is completely described by the symmetric 
adjacency matrix 



A, 



A 



(2) 



All quantitative information about the network geom- 
etry follows from structural properties of the adjacency 
matrix, which can be characterised by a rich variety of 
different measures [9 10II 11I121 . We emphasise that al- 
though we are considering simple graphs, the resulting 
graph topology may be highly non-trivial, which justi- 
fies the term complex network as a contrast to regular 
chains or lattices. 

Recently, it has been suggested that RPs can be al- 
ternatively viewed as a complex network 11311411511611 7 
which captures the geometric skeleton of the attractor 
in phase space (i.e., temporal recurrence relationships 
correspond to spatial proximity relationships). Specif- 
ically, neighbouring observations in phase space are 
represented by mutually linked vertices of a complex 
network, i.e.. 



Aij — Aij(^s) — Rij{e) — 5ij, 



(3) 



where 5ij is the Kronecker delta used here to avoid 
self-loops in the network. We note that this setting cor- 
responds to a specific choice of the Theiler window in 
RQA. As a specific type of networks, the e-recurrence 
networks described by the adjacency matrix Aij are ge- 
ometric graphs [19 20| aka spatial networks [21|, i.e., 
graphs whose vertices are objects in some metric space 
(specifically, the phase space of a dynamical system 
or its reconstruction, e.g., based on time-delay embed- 
ding). More specifically, they can be classified as prox- 
imity graphs [22 1 in arbitrary spatial dimensions. As a 
consequence, given a proper sampling from the consid- 
ered (dissipative) system, e-recurrence networks ap- 
proximate the underlying continuous system and en- 
code relevant geometric information about the corre- 
sponding attractor Moreover, the quantitative analy- 
sis of e-recurrence networks provides complementary 



theoretic quantities are not based on temporal correla- 
tions or explicit time ordering (e.g., line structures in 
the RP) like most classical RQA measures [17 J. 

We note that the idea beyond e-recurrence networks 
is a straightforward generalisation of recent approaches 
in neurosciences [23 24| and climatology [25 26|, where 
the considered systems are approximated by networks 
based on mutual correlations between the individual 
sets of observations measured at certain discrete points 
in physical space. Furthermore, we emphasise that the- 
re are close links to other problems utilising distances 
between objects in some metric space, such as clus- 
ter analysis [22|, dimensionality reduction (e.g., mul- 
tidimensional scaling [27| or isometric feature map- 
ping [281), or set-oriented approaches for identifying 
dynamically invariant objects [29i30j| . 

It should be mentioned that there are multiple ot- 
her approaches to analysing time series by means of 
complex network methods. Specific concepts proposed 
so far include transition networks based on a coarse- 
graining of phase space ll3TI . cycle networks [32|, cor- 
relation networks 1.331 . visibility graphs [34|, and k- 
nearest neighbour ||35l as well as adaptive nearest- 
neighbour networks ['361. The two last methods differ 
from the e-recurrence networks only in the way the lo- 
cal neighbourhood of a vertex is defined. We emphasise 
that the latter class of methods offers the most general 
applicability to a variety of different situations (for a 
detailed discussion of all approaches and their poten- 
tials and limits, we refer to [ 1 7.18] ). The above men- 
tioned methods have already been used for studying 
jgpiplex systems from various perspectives and fields 
or Applications, however, the evaluation of their full po- 
tential is still in the process of exploration. For the 
remainder of this paper, we will exclusively consider 
e-recurrence networks, implying that the results ob- 
tained in this work do not apply to other types of time 
series networks due to their different construction prin- 
ciples. Specifically, even for /c-nearest neighbour net- 
works, the construction principle of which is most sim- 
ilar to that of e-recurrence networks, the local neigh- 
bourhood definition is already so different that our con- 
siderations cannot be easily transferred to this type 
of networks. However, since fc-nearest neighbour net- 
works are based on the original RP definition [2| and 
have some further interesting properties, it will be worth 
considering them in a similar way as done here in fu- 
ture work. 

Recent numerical findings revealed that among ot- 
her measures from graph theory, the local and global 
transitivity properties are particularly well suited for 
identifying and discriminating qualitatively different 
types of dynamics. Specifically, the global network tran- 
sitivity has been demonstrated to provide a good dis- 
criminatory statistics for automatically distinguishing 
between periodic and chaotic dynamics in complex bi- 
furcation scenarios [37|. Moreover, it has been sugges- 
ted that the local clustering coefficient is able to trace 
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II14I17L In this work, we will further elaborate on the 
relationship between attractor properties on the one 
hand, and e-recurrence network properties on the other 
hand, with a special emphasis on dynamically invari- 
ant properties such as fractal dimensions. In partic- 
ular, we will further explore the interrelationships be- 
tween local recurrence rates (i.e., the relative frequency 
of edges a vertex contributes to), the local and global 
transitivity properties of the graph, and the fractal di- 
mension of the underlying attractor. We note that the 
latter results in certain scaling properties of different 
statistical measures that are exclusively based on the 
temporal evolution of the system under study, estab- 
lishing a close link between the dynamics on, and the 
structure of the attractor, which shall be explored here 
from a complex network perspective. We note that sim- 
ilar findings have also been reported independently con- 
cerning links between the structure and function of 
complex networks [111, e.g., in terms of their synchro- 
nisability [38 |, the spreading of failures [39 1, etc. 

This paper is organised as follows: In Sec. [2| we 
review some basic concepts from complex network as 
well as dynamical systems theory needed in this work. 
Some rigorous theoretical results linking the transitiv- 
ity properties of e-recurrence networks with the fractal 
structure of the underlying attractor are provided in 
Sec. [3] and illustrated by numerical findings in Sec. |4] 
Finally, our main results are briefly summarised. 



2 Theoretical background 
2. 1 Network properties 

2.1.1 Direct connectivity properties 

In complex network research, the most important ver- 
tex property is its degree (frequently also referred to 
as degree centrality), i.e., the number of links to other 
vertices in the graph: 



(4) 



Since this measure is extensive, i.e., ki usually increa- 
ses monotonously with increasing N, we prefer to use a 
non-extensive property, the normalised degree or local 
degree density, which can be defined as the probability 
that a given vertex i is linked with any other randomly 
chosen vertex j: 

p, = P{A,, = \). (5) 
For a finite graph, the latter one is estimated by 



Pi 



N 



N -I 



(6) 



We note that for e-recurrence networks, Aij — Aij{e) 
depends explicitly on the recurrence threshold e (i.e.. 



the spatial scale of coarse-graining), so that also fcj and 
Pi are functions of e. In this case, ^(e) corresponds to 
the local recurrence rate RRi{e) of the observed state 
x{ti). In a similar way, we can identify the global edge 
density 



1 



P = 



N{N - 1) 



1 ^ 

I) EA = N^ 



(7) 



of an e-recurrence network with the (global) recurrence 
rate RR observed for the underlying system. Since in 
this case, A^j depends explicitly on e, so do RR as well 
as all graph-theoretical measures derived from the ad- 
jacency matrix. 



2. 1 .2 Transitivity properties 

In general, the term transitivity refers to the reproduc- 
tion of binary relations between mathematical objects. 
Formally, given a set X of objects, a (directed or undi- 
rected) relation R over X is called transitive iff when- 
ever A e X is related to B e X and B is related to 
C e X, then A is also related to C. 

In dynamical systems theory, there are several no- 
tions of transitivity, e.g., metric and topological transi- 
tivity, which describe properties of certain (continuous) 
topological transformations. In contrast, in a complex 
network, the term transitivity is related to fundamen- 
tal algebraic relationships between triples of discrete 
objects [40 41 1. Specifically, in graph-theoretical terms, 
we identify the set X with the set of vertices V, and 
the relation R with the mutual adjacency of pairs of 
vertices. Hence, for a given vertex i £ V, transitivity 
refers to the fact that for two other vertices j,k e V 
with A^j = Ajfc = 1, Ajk = 1 also holds (here as well 
as in the following general considerations, we omit the 
e-dependence of the adjacency matrix for e-recurrence 
networks). In a general network, this is typically not 
the case for all vertices. Consequently, characterising 
the degree of transitivity (or, alternatively, the rela- 
tive frequency of closed 3-loops, which are commonly 
referred to as triangles) with respect to some individ- 
ual vertex or the whole network provides important 
information on the structural graph properties, which 
may be related to important general features of the un- 
derlying system. 

Following the above general considerations, the lo- 
cal transitivity characteristics of a complex network 
are quantified by the local clustering coefficient [42 1, 
which measures the probability that two randomly cho- 
sen neighbours of a given vertex i are mutually linked, 
i.e., 



C, = P{Ajk = l\A.,, = l,Ak = l) 

^ P{A^j = l,A^k = l,Ajk ^ 1) 

P(A„- = l,A,fe = 1) 



(8) 



For finite graphs, the corresponding probability is typ- 
ically estimated in terms of the relative frequency of 
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links between the neighbouring vertices of a given node 
i, i.e., 



C, = 



number of triangles including vertex i 
number of triples centred on vertex i 



(9) 



In contrast, the corresponding correlations between Ci 
and ki, which have been numerically studied in a pre- 
vious paper [17| for different low-dimensional dynami- 
cal systems, have found to be system-specific and often 
not significant. However, there are examples such as 
the logistic map 



ki {ki 



I) 



(12) 



where triple refers to a pair (j, k) of vertices that are 
both linked with i, but not necessarily mutually linked. 

Local transitivity properties translate to global net- 
work properties by sophisticated averaging. In this con- 



text, two different measures maybe distinguished [10 11 
On the one hand, the global clustering coefficient in- 
troduced by Watts and Strogatz Il42l is defined as the 
arithmetic mean of the local clustering coefficients ta- 
ken over all vertices of the network: 



1 ^ 



1 ^ 



(10) 



A potential disadvantage of this measure is that it gives 
equal weights also to vertices with sparse connectivity, 
which can by definition only contribute to few trian- 
gles in the network. On the other hand, the definition 
of the clustering coefficient according to Barrat and 
Weigt [43 44|, which has been later suggested to be 
termed network transitivity [111 , gives equal weight to 
all triangles in the network: 



r 



3 X number of triangles in the network 
number of linked triples of vertices 

X^i j-,fe=l -^jk-AijAik 
J2i,j,k = l ^ij^tki^ — Sjk) 



(11) 



Note that differences between both measures are gene- 
ric and typically persist even for large networks HO J . 



2.1.3 Relationsl^ips between different measures 

As shown above, estimates of both local degree den- 
sity and local clustering coefficient can be written in 
terms of (joint) probabilities for the existence of edges 
in certain parts of a complex network. In terms of e- 
recurrence networks (or, even more general, spatially 
embedded graphs), these different parts may be in- 
terpreted as certain regions in (phase) space. Equa- 
tion ([9| suggests a possibly nontrivial relationship be- 
tween the two aforementioned local network proper- 
ties. For different models of scale-free networks, it has 
be en sh own that ^ k^^ at least for large degrees 
ki II45I46 I. Similar observations have been made for dif- 
ferent real-world networks [47 48 49 1. 

Unlike for the aforementioned examples, for e-re- 
currence networks as spatially embedded graphs, the- 
re is no obvious simple dependence of C, on the den- 
sity of vertices (and, hence, the local degree density). 



with X e 5 c [0, 1] and a e [0,4], where local maxima 
of the degree centrality roughly coincide with maxima 
of the local clustering coefficient (see Fig. jH this fea- 
'I'^Tfe will be further discussed in Sec. |4.1| of this pa- 
per). In general, the local transitivity properties of e- 
recurrence networks do not simply capture density ef- 
fects, but have a distinct meaning in terms of attractor 
geometry, which will be further highlighted in Sec. |3] 

It should be noted that for e-recurrence networks, 
degree centrality, local degree density, and edge den- 
sity are monotonously increasing functions of e by def- 
inition. A similar observation has been made for the 
global clustering coefficient, which follows mainly from 
the fact that parts of the phase space outside the at- 
tractor boundaries get an increasing weight as the prox- 
imity threshold e becomes larger [17|. Note that there 
is no similar general dependence for the local cluster- 
ing coefficient. We will also further address this point 
within the course of this paper. 



2.2 Attractor properties 

The fundamental structural properties of attractors in 
dissipative dynamical systems follow from their invari- 
ant density p{x) associated with the natural measure 
fj. as dfi{x) — p{x)dx. The study of the latter is particu- 
larly interesting for chaotic systems, where the attrac- 
tor has a complex shape in phase space, and allows 
the definition and investigation of dynamically invari- 
ant properties such as entropies, fractal dimensions, 
and related measures of complexity. In this work, we 
are particularly interested in the concept of fractal di- 
mensions and their relationship with the properties of 
c-recurrence networks. 



2.2.1 Global attractor dimensions 

The most common definition of a fractal dimension ta- 
kes a self-similarity property of chaotic attractors into 
account: Given a partitioning of phase space into a set 
of fixed m-dimensional hypercubes of length one de- 
termines the number n{l) of such cubes necessary to 
fully cover the attractor Typically, one observes n{l) ^ 
with some scaling exponent Dq, which is referred 
to as the box-counting or capacity dimension 



Dq = lim 



logn(/) 



f-^a log j 



(13) 
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Fig. 1. Colour-coded representations of estimates of (A) local degree density pi and (B) local clustering coefficient 
C, for the e-recurrence networks obtained from trajectories of the logistic map for different control parameters 
a (N — 10, 000, no embedding, maximum norm, e = O.OScr with a being the empirical standard deviation of 
the respectively considered realisations). (C) Scatter plot between pi and Cj in a; g [0.1, 0.9] (this choice reduces 
the effects of the attractor boundaries) for a ~ 3.9, yielding a rank-order correlation coefficient ps ~ 0.25. (D) 
Magnification of the profile of pi (black) and Ci (grey) in the vicinity of a supertrack function (vertical line) for 
a = 3.9. 



Although Dq is often called "the" fractal dimension, 
there is a multiplicity of other definitions of fractal 
dimensions. These generalised dimensions are related 
with the order-g Renyi entropies 



1 



n{l) 



(geM+U{0}), 



as Il50l 



lin.M) 

i-^o log J 



(14) 



(15) 



with Dq^ < Dq^ for qi > q2. Special cases of this defini- 
tion include Dq, the information dimension Di (where 
5*1 corresponds to the Shannon entropy), and the cor- 
relation dimension D2. The latter can be equivalently 



defined using the correlation integral IISDl 



C{e)= / dp{x) / dp{y)e{e~\\x^y\\), 



as 



Dy = lim 



logC(£) 
loge 



(16) 



(17) 



Given only a finite set of sampled points on a specific 
trajectory (i.e., a time series) of the system, an unbi- 
ased estimator of C (e) is given by the correlation sum 



C(e) = lim 



1 



N{N 



N 



e{s-\\x,-xj\\), (18) 
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which corresponds to the global recurrence rate RR 
(and, hence, the edge density p of the associated e- 
recurrence network). If the supposed power-law beha- 
viour C(e) cx of the correlation integral applies, one 
may usually find that also p{e) cx e^^ for s e [smm, £max], 
such that D2 can be estimated as 

^^^dl0^p(£) (19) 

aloge 

The finite size of this scaling interval is caused by the 
fact that (i) for small e, the number of recurrences (i.e., 
neighbours in phase space with respect to the consid- 
ered e-distance) is too low to observe the correct slope 
of log p = /(loge) with sufficient statistical confidence 
(finite-resolution limit of a finite time series), whereas 
(ii) for large e, parts of phase space that are not covered 
by the attractor (and therefore do not contain any ob- 
servations) are successively included in the considered 
e-balls, which leads to a saturation of p. Moreover, in 
the latter case, more and more redundant information 
is included in the estimation (e.g., due to the consid- 
erations of points that are nieghbours in phase space 
just because they are also close in time). Therefore, in 
most available methods for estimating D2 from time 
series, the consideration of a sufficiently large ensem- 
ble of different values of e is necessary to obtain proper 
estimates. 

A practical alternative is considering scale-local di- 
mensions Dq{s) II52I . i.e., numerical estimates of Dq 
obtained from the respective relationships evaluated 
for only one fixed value of e. Note, however, that these 
quantities often do not provide proper estimates of the 
true dimension. 



2.2.2 Pointwise dimensions 

The defining equation ( |17[ l of the correlation dimen- 
sion has been formulatedglobally as a property of the 
whole attractor. However, the correlation sum is de- 
fined as an unweighted average of contributions from 
all observed state vectors x^, i.e., 

1 ^ 

Cie) ^ ^Y.^^{e) (20) 

1=1 

with 

1 N 

= E 0{e-\\x,-x,\\)=p.{e) (21) 

coinciding with the local recurrence rate (local degree 
density pi) of Xi (the same considerations apply to the 
correlation integral C(e) in Eq. ([l6)l). Motivated by this, 
one defines the pointwise dimension ||53l 

i^,(x) = lim ^°g^^^-("» , (22) 
' e^O loge 



where p.{B^{x)) is the measure of a ball of radius e cen- 
tred at X. Commonly, this pointwise dimension is esti- 
mated as 

Dp{xi) = for e e [e™„,£max]- (23) 

a log e 

If Dp{x) exists (note that convergence cannot be ex- 
pected in general), it is independent of x for almost all 
x with respect to the invariant measure p. Since un- 
like for C(e), no two-point correlations are involved in 
the local recurrence rate, the pointwise dimension pro- 
vides a local estimate of the information dimension Di 
rather than D2. As for the (global) correlation dimen- 
sion, the practical estimation of pointwise dimensions 
may be challenging since the limit e — > can hardly be 
assessed with a limited number of data available. 



2.2.3 Lyapunov dimension 

The Lyapunov (or Kaplan- Yorke) dimension is based 
on the temporal stretching and folding characteristics 
of the dynamical system under study. Let {Ai, . . . , A™} 
denote the spectrum of Lyapunov exponents of the sys- 
tem ordered such that Ai > A2 > • • • > Xm, and n be the 
largest integer such that 

n 

Xj > 0. (24) 

Then the Lyapunov dimension or Kaplan-Yorke dimen- 
sion Dl is defined as 1,53.54 55) 

1 " 

^^ = "+U lE^J- (25) 

The Kaplan-Yorke conjecture states that Di = Dl for 
typical attractors. 

The Lyapunov dimension is typically considered a 
global measure as it applies to the whole attractor ex- 
cept for a set of measure zero. In contrast, the local 
Lyapunov dimension Dl{x) can vary across the attrac- 
tor [56 57 1. It has (so far) mainly been studied for maps 
F and can be obtained from Eqs. ( |24[ |25| l after the re- 
placements Dl — > Dl{F,x) and A^ — >■ ln(aj(x)), where 
Gj (x) are the ordered singular values of the map's Jaco- 
bian DF{x). The most important property of Dl{F,x) 
is that it provides an upper bound on the box-counting 
dimension Dq, i.e., 

< maxDi(f,a;). 

X 

The local Lyapunov dimension Dl{F^,x) of the k-ih. it- 
erate of the map F converges to the Lyapunov dimen- 
sion Dl as fc 00 [56|. 
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2.3 Links between network and attractor properties with x e S and e > are p-measurable, where 



We have already argued that (given a proper sampling 
of the considered trajectory) the local degree density of 
an e-recurrence network (i.e., the local recurrence rate) 
allows defining an estimator of the invariant density 
p{x) of the underlying attractor. The latter is known 
to be fundamental for both the structure of, and the 
dynamics on the attractor, whereas a similar state- 
ment applies to the connectivity of vertices and the 
overall structure of a network. These established links 
suggest that there might well be more possible inter- 
relationships between dynamical system and network 
properties. This is evident for the estimation of the cor- 
relation dimension from the scaling of the edge den- 
sity with varying s. In addition, as it will be shown 
in a complementary paper Il58l . there are cases where 
the degree distribution of an e-recurrence network in- 
cludes a scale-free part, the characteristic exponent of 
which can be related to the pointwise dimension of 
the underlying system. We note that similar observa- 
tions have been made for another type of complex net- 
works generated from time series, the so-called visibil- 
ity graphs, which also show a power-law decay of the 
degree distribution for certain fractal time series, the 
exponen t of wh ich is directly related to the Hurst pa- 
rameter II59I60L 

Another important relationship between the geo- 
metric properties of a dynamical system and the tran- 
sitivity properties of their induced e-recurrence net- 
works can be inferred from the theory of random geo- 
metric graphs Il6l]| . Specifically, it has been shown that 
for such graphs, the expected global clustering coeffi- 
cient depends on the (integer) dimension of the metric 
space in which the considered graph is embedded. In 
the remainder of this paper, we will generalise this re- 
sult to arbitrary non-integer spatial dimensions and 
discuss the resulting implications in some detail. 



3 Continuous clustering coefficient and 
clustering dimension 

In the following, we develop and successively apply a 
general theory linking the local as well as global tran- 
sitivity properties of e-recurrence networks with geo- 
metric attractor properties. For this purpose, we de- 
fine a novel notion of fractal dimension based on these 
transitivity properties, and compare these new defini- 
tions with several other existing measures described 
in the previous section. 

3.1 Continuous measures for transitivity 

3.1.1 General theory 

Let p he a probability density on some set S c R™, 
such that all closed e-balls 

B,{x)^{yeS:d{x,y)<e} 



d{x,y) 



max \x., - 

i—l,...,m 



is the maximum metric and S is the topological closure 
of S. Then we define the continuous e-degree density 
and continuous e-clustering coefficient of any point x £ 
5 as 



p(x;e) 



C{x;e) 



P{y) dy= dfj.{x) 

B^x) JB,(x) 

IIb,{x) P(y)p(^)'9(£ - d{y, z)) dy dz 
p{x]eY 



(26) 



, (27) 



the latter being the probability that two points y and z 
randomly drawn according to p are closer than e given 
they are both closer than e to x. As a global measure, 
we define the continuous e-transitivity of S as 



Tie) = 



p{x)p{y)p{z)e{e - d{x,y)) x 



x6'(e — d{y, z))0{e — d{z, x j) dxdydz 



p{x)p{y)p{z)0{e - d{x, y)) x 



x6'(e — d{z,x)) dx dy dz 



(28) 



which is the probability that among three points x, y, z 
drawn randomly according to p, y and z are closer than 
e given they are both closer than e to x. 



3. 1 .2 Examples: Dynamical systems defined on the 
unit interval 

If S is the unit box [0, 1]™ and p is the uniform den- 
sity p{x) = 1, one can easily see that C{x:e) = T{e) = 
(3/4)™ for all points a; e [e, 1 — e]™. This is because for 
m = 1, a randomly chosen y e B^{x) — [x — e,x + e] 
has on average three quarters of within its own 

B^iy), namely half of B^{x) for y — x±e and the whole 
B^{x) for y — X. For m> 1, the integral in the numera- 
tor of Eq. (|27]l decomposes as 



C(.x;e) 



1 



0{e-\7j-C\)dr,dC 

i^l — x.-e (29) 

2\m //o-^2m 



n 



(2e^2m 

(3e")™/(2e)"™ = (3/4) 



The same is true whenever x is in the topological 
interior of S, p is sufficiently smooth in a neighbour- 
hood of X, and e is small, because then the situation 
looks locally approximately the same as for the uni- 
form density on the unit box. For example, if S and p 
are the main attractor and invariant density of the lo- 
gistic map ([12| at a = 4, for which S = (0, 1) and El 



p{x) 



1 



TT^/x(T 



(30) 
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then one can also show that for small e, 

2e 



and 



so that 



p{y)p{z)e{s - d{y, z)) dy dz 



C{x;e) 3/4 



as £ -> 0. For x e {0,1}, we get C{x;e) = 1 for all e, 
since then all pairs of points in B,. (x) are on the same 
side of X and thus have a distance < e. Note that these 
results are consistent with recent findings described in 
|[nj(see Sec.liUl. 

In general, the stronger a smooth p varies inside 
B^{x), the more C{x; e) exceeds its lower bound (3/4)™. 
More precisely, the larger m, the closer C(a;; e) is related 
to the 2nd-order Renyi entropy of p restricted to B^{x), 
since the Heaviside function then acts more and more 
like a delta function and the denominator of C{x;e) 
becomes approximately proportional to Jg (^^^(y)^ dy. 
The phenomenon that the clustering coefficient is in- 
fluenced by the local density variability rather than by 
the local density level was also observed in a different 
context for climate networks [63 1. 

We have to emphasise that the above result C(a:;; e) ~ 
T{e) — (3/4)™ is only valid when using the maximum 
norm for defining distances between points in phase 
space. For other choices of the metric d, e. g., the Eu- 
clidean one, there will in general not be a similarly 
simple exponential relationship, although it is some- 
times at least asymptotically exponential. For exam- 
ple, for the uniform density on the unit box and using 
the Euclidean metric, we find 



C{x;e) = 1 - 



mr{m/2) 



2^r{{m + l)/2) _ 



2F1 



1 1-m 3 1 



2' 



1 /I — m TO + 1 m + 3 1 

2F1 



m + 1 



2 ' 2' 4 
2 '4 



TTTO V 4 



(m+l)/2 



0.862™/3, 



(31) 



where 2^'i( ) is the hypergeometric function and the 
exponential fit is only good for m > 40. The latter result 
can be obtained from the expressions for the volumes 
of m-dimensional spheres and their intersections and 
is consistent with previous findings in (HT). 



3.2 Clustering dimensions 

3.2.1 General theory 

The aforementioned observations motivate the defini- 
tion of two new local and two new global measures of 



dimension for general S and p, namely the upper and 
lower clustering dimension of S at x, 



Dq{x) ~ limsup 



logC(a;;£) 



0" log(3/4) 



and Dq{x) = liminf ■ 



>o log(3/4) ' 
and the upper and lower transitivity dimension of S, 

logr(e) 



(32) 
(33) 



DVj- = limsup 



and Dq- = lim inf 



0" log(3/4) 
. . , logr(e) 



log(3/4)- 



(34) 
(35) 



Note that [x] and D\. {x) need not be continuous 
in x. For example, in the logistic map with a — 4, the 
border points x e {0, 1} have C{x; e) = 1 (see above), so 
that D'c{x) = D\,{x) jumps from 1 to for x — > 0, 1. For 
a < 4, such discontinuities also occur in the interior of 
S, because inside Bs{x) for a point x on a supertrack 
function [64 65 1 of the logistic map, p has asymptoti- 
cally a power-law form 

p{y) ^ \y-xy 

on one side of x, and is approximately constant on the 
other side of x [55 1 . To see this, assume that the k- 
th iterate of the map, /''', has a local maximum at x, 
f^{x) = y. Then, all points close to x get mapped by 
to points rather close to y but smaller than y. This 
explains why the density then has a peak at y, which 
diverges only to the left side of y. The power-law de- 
cay can be seen from the local quadratic approxima- 
tion of at X, since that has a slope linear in the 
distance from x, and the density transforms using the 
inverse of the slope, i.e., l/( distance from x). As a con- 
sequence, for e — > 0, almost all mass inside Bs{x) is on 
the power-law side and, hence, almost all y,z ^ B^{x) 
have distance < e, yielding C{x;e) — > 1 for e 0. 
Hence, D^{x) = D\,{;x) = at all points on the su- 
pertrack functions of the logistic map. This behaviour 
can be seen in Fig. [l] where the supertrack functions 
are clearly visible in terms of pronounced maxima of 
C{x) that have been numerically estimated from the 
e-recurrence networks of sample trajectories 1,66 J . 



3.2.2 Examples: Self-similar sets 

If S is highly self-similar, D'c{x) and D^^{x) can also 
differ considerably, as can be seen from the example in 
which p(a:;) = 1 and S* = {x e [-1, 1] : \x\ e Sc}, where 
Sc C [0, 1] is the Cantor set (= numbers that have a 
ternary expansion in which all digits are either or 
2). In that case, it is easy to see from the symmetry 
that C(0; e) oscillates between 3/4 and 1 for e — ^ 0, with 
C(0;e) = 3/4 if e = 1/3" for some integer n > 0, and 
C(0;e) 1 if e 2/3" for some integer n > 1. Hence, 
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D^iO) = 1 and £)^(0) = 0. Similarly, for almost all x e 
S, one can show that C{x; e) is infinitely often 7/8 and 
1 for e ^ 0, so that D'^{x) > log(7/8)/log(3/4) « 0.464 
and D\.{x) = 0. Only for the countably many points 
of the form x — ±a/3" with integers a,n > 0, we get 
Ce{x) = 1 for all e < 1/3" and, hence, D^{x) = D\.{x) = 
0. In the same example, also Dj- and Dj- differ since 
T{e) = 1 whenever e = 1/3" with integer n > and 
T{e) = 11/13 whenever e = 5/3" with integer n > 1, 
with intermediate values for other values of e, so that 
= log(ll/13)/log(3/4) w 0.581 and = 0. 

The same values are obtained when 5' is a version 
of the Cantor set in which, starting with the unit in- 
terval, iteratively the middle fraction of relative width 
1 — 2q! is removed from every remaining interval, for 
a < 1/3, where a 1/3 gives the standard Cantor set. 
Other measures of dimension, however, have values 
depending on a, e.g., the box-counting, information, 
and correlation dimension of S and the pointwise di- 
mension of almost all x £ S are all — log 2 / log a, which 
varies between log 2/ log 3 « 0.631 and 0. In particu- 
lar, this shows that both (x) and D}^ can be either 
smaller or larger than all those classical measures of 
dimension. This finding is supported by further analyt- 
ical as well as numerical results (see, e.g.. Sec. 4.2. 2[ ). 

As a further example, let us consider the generalised 
baker's map 



Xn+l 
Vn+l 



(1 - Ab) + Xbx„ 



< a, 
> a, 



Vn/a, Vn < a, 

iVn ~ a)/(l - a), Vn > a 



(36) 



with a < 1, Aa, Ab > and Aq + Ab < 1 II53I55I . which 
yields a transformation of the unit square [0, 1] x [0, 1]. 
If S is the attractor of the symmetric version of this 
map with a = 1/2 and Aq = Ab = 1/4), it is a carte- 
sian product of a Cantor set similar to those discussed 
above and the unit interval. Hence, for almost all x e 
5*, C{x;e) attains 7/8 • 3/4 and 1 • 3/4 infinitely often 
for £ ^ 0, so that D'^{x) > log(7/8)/log(3/4) + 1 « 
1.464 and D'^{x) = 1. Similarly, we obtain that = 
log(ll/13)/log(3/4) + 1 « 1.581 and Dij- = 1. 

The above examples suggest that it might be inter- 
esting to take the difference — Dij- as a measure 
of self-similarity of an attractor. As a counter-example, 
one may consider again the chaotic attractor of the lo- 
gistic map at a = 4, for which D^^'\x) = Z?^'' = Dr ^ I 
is consistent with Dp{x) ~ Di — m — I as one would 
expect for this non-fractal case. 



4/7 



5/7 



Fig. 2. Intermediate step S2 (black squares) of the 
fractal construction to show that can exceed m. 



2.24 > m since C{0;e) = 21/40 < 3/4 for all e e (0, 1). 
This is because of p{0) = 0, so p{y) cannot be considered 
approximately constant on B^{0) for small e, although 
p is smooth. The best upper bound for (x) in terms of 
m is m log 2/ (log 4 — log 3) w 2.41 m. To see this, we split 
B^{x) into 2™ generalised "quadrants" corresponding 
to the possible vectors of signs yi — Xi, i = 1 . . . m. 
Note that all pairs y,z e B^{x) that are in the same 
generalised quadrant have d{y, z) < e, and the prob- 
ability of y, z being in the same generalised quadrant 
is at least 1/2™ no matter how p{B^(x)) is distributed 
over these regions, so that C{x\£) > 1/2™. That this 
bound is sharp even for arbitrarily smooth p can be 
seen from the set of examples with S = [—1,1]™ and 
Pi^) = ni"i(2« + l)2;f°/2 with integers a, m > 0, for 
which C(0; e) 1/2™ as a ^ 00. 

That might also exceed m can be seen from an- 
other Cantor-like example. Starting with the unit cube 
So = [0, 1]™, replace this cube by 1+2™ smaller cubes of 
the form n,=i[3/7 + s„ 4/7 + Si] with si = • • • = s,„ = 
or Si e {—3/7, 3/7} for all i, i. e., one small cube located 
at the center of the original cube, and 2™ cubes fit into 
the corners of the original cube, giving a set 5*1 c Sq. 
Repeating the same replacement infinitely often with 
each cube gives a descending sequence of sets 5*^. Fig- 
ure [2] shows the result 5*2 of two iterations for m = 2. 
The intersection S = Pli^i = ^oa is a Cantor-like 
fractal for which T{e) = (7 • 2™ + l)/(4™ + 6 • 2™ + 1) < 
(3/4)™ whenever m > 5 and e = 4/7" for some inte- 
ger n > 0, hence D^{S) > m. We do not know whether 
even Dl^{S) can exceed m but conjecture that this is 
impossible (this conjecture is supported by further nu- 
merical results in Sec.|4]l. 



3.2.3 Clustering and topological dimension of 
phase space 

In some cases, D^{x) and D\,{x) may even exceed the 
non-fractal dimension m of the surrounding space. For 
example, for S = [—1,1] and p{x) — 3x^/2, Dj.{x) — 
D^{x) = m = 1 for all x ^ Q, but D\.{Q) ^ D^{0) w 



3.2.4 Local clustering and Lyapunov dimensions 

For a homogeneous fractal the pointwise dimension as 
the most traditional other local dimension measure by 
definition equals the global measure of information di- 
mension for all X except for a set of measure zero, i.e.. 
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Dp{x) = Dp = Di II53II . However, many chaotic attrac- 
tors have a complex internal structure (e.g., embedded 
objects of measure zero with deviating pointwise di- 
mensions), such as the Rossler [67 1 and Lorenz sys- 
tems [67 68|. In such cases, D^{x) and Dc{ x) ca n show 
considerable regional differences (see Sec. 4.3 for the 
Rossler system). Beside the pointwise dimensions, the 
only other local measure of dimension with this prop- 
erty that we are aware of is the local Lyapunov dimen- 
sion Dl{x) fb^, which is based on the local contraction 
of a map. Thus, it might be worthwile to compare these 
measures. 

For the logistic map, e.g., Dl{x) = lior x ^ [3/8, 5/8] 
and Dl{x) = otherwise, while D^{x) = D^(,{x) = 1 for 
X e (0, 1) and D^{x) = D\,{x) = for x e {0, 1}. For the 
generalised baker's map, we will provide some results 
in SeciXll 



3.3 Estimation of ciustering dimensions 

In general, if S is an attractor of an ergodic dynami- 
cal system with invariant density p, we can estimate 
C{x{ti);£) and T{£) from a finite sample x{ti) G 5 of a 
trajectory in S sampled at time points ti,...,t]s, using 
the standard (sample) clustering coefficients (Eq. (|9]l) 
and transitivity (Eq. (111). 

By the weak law of large numbers, Ci (e) and T(e) 
converge in probability to C(.T(ij);e) and T(e) as iV — > 
CO, for each e > and a general choice of time points 
ti (e. g., regularly spaced with a time step that is co- 
prime with all periodic orbits' periods). In other words, 
Ci (e) and T(£) are statistically consistent estimators of 
C{x{ti);e) and T{e). Our four new dimension measures 
can then be estimated from a sufficiently long sampled 
trajectory as 



^c,i — 



and b 



T 



logC,(e) 

max 7— T7T: 

eee log(3/4)' 

. logC»(g) 
™G£ log(3/4)' 

iogr(£) 

max 7— T7T: 

eef log(3/4)' 

. logr(£) 

mm — 

eG£ log(3/4) 



(37) 



(38) 



(39) 



(40) 



for a suitably large set E of different values of e that 
are as small as possible while still providing for suf- 
ficiently large values of hi(e). We emphasise that the 
above set of equations is only feasible if Ci^e) > and 
T(e) > 0, respectively. For all other points in phase 
space represented by a vertex i, the aforementioned di- 
mension measures cannot be defined in a meaningful 
way. 

Note that h{e) has a binomial distribution and be- 
comes (for e — > and large N) asymptotically Poisso- 
nian with mean and variance 11211 

\ ^ Np{x)e°'-^''\ 



where Dp again denotes the standard pointwise dimen- 
sion. For a given fci(e), Ci{e) has asymptotically the 
mean Ce{x{ti)) and a variance proportional to fc,(e)~^, 
so that also Ci(e) has asymptotically a variance pro- 
portional to ki{e)^^. For this reason vertices with low 
degree are unlikely to yield reliable estimates of local 
clustering dimensions (see below). 

The quotient max£/min£ should exceed the mag- 
nification factor of any suspected self-similarity (i.e., 
3 in the Cantor set), see Fig. |9p for an example. In 
general, it has been established that there is no gen- 
erally applicable rule for choosing e for computing re- 
currence network properties [16|. In fact, the values 
of e that may provide feasible results are restricted 
by practical considerations, i.e., the available sample 
size determines the smallest possible spatial dimen- 
sions of structures to be resolved by network-theoretic 
measures, which directly relates to the minimally pos- 
sible e, whereas large e do not allow obtaining informa- 
tion on the geometric fine structures of the attractor. In 
typical situations, a reasonable trade-off can be found 
for edge densities below 5%. With respect to the joint 
estimation of upper and lower clustering/transitivity 
dimensions, our numerical studies (see Sec. |4]) reveal 
that for attractors with a self-similar structure, these 
measures typically alternate between lower and upper 
bounds as e is varied (cf Fig. [9j3). According to this 
observation, we suggest as a rule of thumb that the 
range of e should be chosen so that both upper and 
lower limits can be identified from at least two distinct 
intervals of e, respectively. Note that in contrast to the 
estimation of other notions of dimension (like with the 
Grassberger-Procaccia algorithm in the case of corre- 
lation dimension [69 1), it is not necessary to estimate 
the slope from a double-logarithmic plot here. 



4 Examples 

In the following, we illustrate our previous consider- 
ations by numerical as well as further analytical re- 
sults obtained for some benchmark examples of both 
low-dimensional maps and time-continuous dynamical 
systems. 



4.1 Logistic map 



4.1.1 The a = 4 case 



For the logistic map ( 12 1 at a = 4, the knowledge of the 



invariant density p{x) of the main attractor (Eq. (301) 



allows deriving analytical expressions for quantities 
such as local degree density and local clustering co- 
efficients [17|. In the latter case, one has to consider 
the simplification that the denominator in Eq. (|8]l fac- 
torises, i.e., the probabilities of two randomly chosen 
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Fig. 3. Analj^ical values (upper panels) and numerical estimates (lower panels) of the local degree density 
p{x;e) (left) and the local clustering coefficient C{x;£) (right) for the logistic map at a = 4. Numerical results 
have been obtained for one realisation of the system with N = 10, 000 points. Diagonal lines indicate regions 
that are affected by information taken from outside the attractor. 



vertices j and k to be in the e-ball around Xi are inde- 
pendent and equal: 



P{A,, = l,A,k = 1) ^ P{A,j = l)P{A,k = 1) 
^P{A,,f^pl 



(41) 



Since the chaotic attractor is bound to the inter- 
val (0, 1), a careful treatment of the resulting integra- 
tion boundaries reveals that spatial variations in all 
measures can be understood as being originated in the 
coverage of phase space regions outside the attractor 
(i.e., outside the interval (0,1)) rather than being a di- 
rect effect of the local degree density The correspond- 
ing results obtained in [ 17 1 are briefly summarised in 
Tab. [l] The analytical considerations explain the nu- 
merical findings concerning the e-dependence of the 
global clustering coefficient C as well as the x-depen- 
dence of the local clustering coefficient C, for fixed e in 
an excellent manner [17|. Figure [3] demonstrates that 
this actually holds for both local degree density and 
clustering coefficient and for all x and e. 

From the dependence of the invariant density and, 
hence, the local degree density on both x and e, one 
may conclude that the pointwise scale-local dimension 
shows a similar dependence. However, when looking 



in more detail at estimates of the actual pointwise di- 
mension Dp{x) obtained from the scaling exponent of 
the local degree density (Fig. one finds that as ex- 
pected, both Dp{x) and the scale-local (fixed e) esti- 
mate Dca{^) = logCi(e)/log(3/4) approach values close 
to 1 in a broad range in the middle of the chaotic at- 
tractor, i.e., in some interval that is not influenced by 
the attractor boundaries. The smaller e and the larger 
N, the better the pointwise convergence of Dp{x) — > 1 
and Dc.i{^) ^ 1 obtained by numerical calculations for 
this region. We emphasise that the profile of the point- 
wise dimension is much smoother than that of the local 
clustering dimension, which follows from the fact that 
a variety of different values of e has been used in the 
estimation of Dp{x), while only one e had to be consid- 
ered for Dc,i{s), resulting in a larger variance of the 
estimate. 

Close to X — e and 1 - e, the corresponding esti- 
mates become larger than the theoretical upper limit 
of 1, whereas for x — > 0, 1, Dp{x) — ^ and Dc.iis) ~^ 
as expected (see Sec. |3.2.l"| . The observed overshooting 
of the estimates Dp{x) close to x = £,l — e results from 
a loss of convergence of the estimator, which is under- 
lined by Fig. kp. Moreover, the scatter plot of Ci versus 
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A 1 




C 1 



Fig. 4. Properties of the e-recurrence networks ob- 
tained for one realisation of the logistic map at a = 4: 
Point estimates of (A) Ci(e) and (B) Dp{xi) as well 
as (D) the goodness-of-fit of the linear regression 
of log/Si vs. loge used for the estimation of Dp{xi) 
(Eq. pS])) in dependence on x. In addition, the rela- 



tionship between estimates of the local clustering co- 
efficient Ci and pointwise dimension Dp{xi) is shown 
for different choices of the edge density p in panel (C). 
For the estimation of network properties, N = 10, 000 
points have been used with an edge density of p ^ 0.01, 
while estimates of the pointwise dimension have been 
obtained from N — 10^ points. The initial 100 points 
have been removed to avoid transient behaviour. 



Dp{xi) for different e (Fig.Uj^) demonstrates that apart 
from the regions close to tne attractor boundaries, all 
points with an e-neighbourhood that lies completely 
within (0, 1) are characterised by Dp{xi) w Dc,i{^) ~ 1 
as expected. For the regions suffering from boundary 
effects (which are particularly pronounced in Fig. |4p 
due to rather large choices of e, cf the discussion in [ 16 1), 
there is still a clear (but nonlinear) dependence be- 
tween both measures. Note that the two-band struc- 
ture in the scatter plot between both measures results 
from some numerical effects due to a slight asymmetry 
between the density close to the two attractor bound- 
aries, which is expected to vanish for higher N. 



4.1 .2 Bifurcation scenario 

It has already been shown that for the logistic map, re- 
gions with a high invariant density, which are typical 
for the supertrack functions, coincide with local max- 
ima of both local degree density p{x) and local cluster- 
ing coefficient C{x) (Fig.[l]). For the local degree density, 
this observation is related to the fact that supertrack 
functions correspond to accumulation points of iterates 
of the map, in the vicinity of which trajectories tend 
to stay for a finite amount of time since they are only 
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Fig. 5. Profile of the invariant density p{x) (blue) and 
estimates of the two network measures pi (black) and 
Ci (red) for the logistic map at a = 3.9 (parameters as 
in Fig.[l]l. The positions of the first 12 supertrack func- 
tions are indicated by vertical lines. Note that unlike 
the invariant density, the degree density is not a prob- 
ability density, i.e., not normalised. 



weakly repulsive in comparison with the usual expo- 
nential separation rate of the map llTOl . A correspond- 
ing reasoning for C{x) has already been discussed in 
Sec. dm 

From the application point of view, one could ask 
whether p{x) or C{x) are better suited for approximat- 
ing the (possibly unknown) location of supertrack func- 
tions in a map. The specific supertrack shown in Fig.jlJ) 
suggests that the local maximum of C{x) approximates 
the theoretical location better than that of p{x). This 



observation is related to the argument from Sec. 3.2.1 
that the invariant density p{x) is approximately con 
stant on one side of a supertrack function of the map, 
but decays like a power-law with increasing distance 
on its other side. As a consequence, estimating p{x) 
with a finite e can be expected to result in a bias of 
the local maximum of p{x). In order to study the gen- 
erality of the latter finding. Fig. [5] shows the complete 
profile of p{x) and C{x) for a = 3.9. One finds that at 
least the most pronounced interior maxima of the lo- 
cal clustering coefficient indeed coincide very well with 
supertrack functions of low order, whereas the corre- 
sponding maxima of p{x) are somewhat shifted from 
the known locations of the supertracks. We emphasise, 
however, that these shifts are directly related to our 
choice of e (the same holds for the differences between 
the estimated Ci and the theoretically predicted value 
of 1 at the supertracks) and vanish in the limit — > oo, 
£ — > 0. For further higher-order supertrack functions, 
the numerical coincidence of maxima of C{x) with the 
supertracks is even weaker, which is also a result of 
the finite sample size and insufficient spatial resolu- 
tion. Using longer realisations and smaller values of e 
(not shown) yields a more reliable profile and, hence, 
improves the skills of both vertex properties for local- 
ising supertrack functions. 
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4.2 Two-dimensional maps 

For the logistic map discussed above, the chaotic at- 
tractors often cover simply connected subintervals of 
(0, 1), with the possible exception of a countable set 
of isolated points on supertrack functions, which has 
measure 0. In contrast to this case, there are numerous 
examples of chaotic maps that have attractors with a 
pronounced fractal structure. Following our theoretical 
considerations from Sec.|3] it is interesting to study the 
behaviour of clustering and transitivity dimensions for 
such maps and compare it with the classical concepts 
of pointwise and local Lyapunov dimensions (see, e.g., 

CD). 




. P,-0-02 



■ 1.0 B 0.4 



1.5 



0.9 
0.8 
0.7 
0.6 
0.5 




0.2 



-0.2 
-0.4 











bp 









-0.5 0.5 



I 



I 



1.5 



2.06 
1.65 
1.24 
0.83 
0.42 



0.5 




4.2.1 Henon map 

As a first example, we consider the Henon map 172 

„2 



(42) 



with the canonical parameters a ~ 1.4 and h = 0.3. The 
chaotic attractor of this map has a fractal structure, 
being smooth in one direction and a Cantor set in the 
other. Numerical estimates of the correlation dimen- 
sion yield D2 = 1.42 ± 0.02 I^. The pointwise dimen- 
sion of the attractor has been extensively discussed in 
the framework of multifractal chaotic attractors and 
unstable periodic orbits ITTIII . 

Figure [6] shows colour-coded representations of the 
local clustering coefficients and pointwise dimensions. 
Unlike for the logistic map, we find no clear relation- 
ship between both measures. A pronounced exception 
are the tips of the attractor, which locally represent 
zero-dimensional structures (Dc{x,y) — > 0, Dp{x,y) — > 
0). The strong differences between the respective mea- 
sures, which can be found in large parts of the attrac- 
tor, seem to be a consequence of the specific filamental 
structure of the attractor. In fact, the Henon attractor 
has a fractal support, which results in rather specific 
topological and metric properties fTSl. 

A similar inconsistency can be observed when com- 
paring the upper and lower clustering dimensions with 
the local Lyapunov dimension (see Fig. [7]l, where no 
clear statistical relationship seems to exist as well. Nev- 
ertheless, the local clustering dimensions behave in 
the expected way: the upper clustering dimension 
is smaller tha n the upper bound 2.41™ « 4.82 dis- 
cussed in Sec. |3.2.3| for almost all vertices, with only 



Fig. 6. Properties of the e-recurrence networks ob- 
tained for one realisation of the Henon map at a = 1.4 
and b = 0.3: Colour-coded representation of point es- 
timates of (A) Cj and (B) Dp{xi) as well as (D) the 
goodness-of-fit for the estimation of Dp{xi) in de- 
pendence on X and y. In addition, the relationship be- 
tween estimates of the local clustering coefficient Ci 
and pointwise dimension Dp{xi) is shown for different 
choices of the edge density p in panel (C). For the esti- 
mation of network properties, N = 10, 000 points have 
been used with an edge density of p — 0.01, while esti- 
mates of the pointwise dimension have been obtained 
from N = 10^ points. The initial 1000 points have been 
removed to avoid transient behaviour. 

specific location, e.g., close to the tips of the attractor 
bands. 

As a consequence of the aforementioned observa- 
tions, very long realisations are typically required to 
numerically capture the local features of the chaotic 
attractor of the Henon map with reasonable confidence. 
This is further underlined by the transitivity dimen- 
sions: the larger N, the better the estimated values of 
this measure obtained for fixed e approach stationary 
values corresponding to the upper and lower transitiv- 
ity dimensions (Fig.jSjl. In contrast, for too small N, we 
observe significant deviations from the asymptotically 
estimated values, which becomes particularly impor- 
tant for £ — ^ 0. 



very few exceptions corresponding to vertices with low 
degree (Fig.^B). A similar observation is made for Dj^, 
which is always smaller than the dimension of the sur- 
rounding space (m — 2) and shows non-zero values for 
the vast majority of vertices (Fig. |7]\). Most vertices 

with zero values of have very low degree as well. Sec. 3.2.2 that for almost all {x,y) e 5*, D^{x,y) 
pointing to a purely statistical effect. However, there 
are some exceptions such as vertices with some very 



4.2.2 Generalised baker's map 

For the symmetric version of the baker's map (Eq.|36|l, 
one finds that for y 7^ 1/2, the Jacobian is given as 
diag(l/4, 2) with singular values ai = 2 > 02 = 1/4. 
Hence, the local Lyapunov dimension is y) = 1 + 

1/2 since ai > 1 > aia2 and 0102^'^ = 1. For the clus- 
teri ng and transitivity dimensions, we have shown in 

> 



1.464, DVj- « 1.581, and Dj.{x,y) = = 1- The latter 
results are confirmed by numerical calculations, the 
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Fig. 7. Relationships between point estimates of (A) 
lower and (B) upper clustering dimensions on the one 
hand, and local Lyapunov dimensions on the other 
hand, for one realisation of the Henon map at a = 1.4 
and b = 0.3 estimated from N — 30, 000 data points 
(initial condition {x,y) = (0,0), the first 1,000 itera- 
tions have been removed from the trajectory to avoid 
transient behaviour) using 100 equally spaced values 
of e in the interval [0.005,0.1]. Red dots correspond to 
vertices with low degree h < 5, for which the obtained 
estimates of D^'l are hardly significant. 



results of which are summarised in Fig. [9] It is notable 
that the estimated values of the transitivity dimension 
roughly coincide with the theoretically predicted upper 
and lower bounds. However, there are examples where 
these bounds are exceeded. We identify two possible 
reasons for such behaviour: too large e or (for small 
e) too small N, i.e., finite-scale and finite sample size 
effects. Due to the resulting outliers obtained when 
varying e, numerical values for the upper (lower) clus- 
tering dimension typically have a positive (negative) 
bias with respect to the theoretically predicted values, 
which is nicely illustrated by Fig. [9]C-F. As a statisti- 
cal estimation effect, this bias is more severe for local 
dimensions, since the variance is larger than for tran- 
sitivity dimensions (see above). However, a bias also 
exists for the transitivity dimensions, where it is just 
smaller (e.g., see the overshooting in Fig. [9^). 




0.04 0.06 

e 



Fig. 8. Estimation of the transitivity dimensions i)^' 
of the Henon map at a = 1.4 and b — 0.3 (one real- 
isation with initial condition {x,y) — (0,0), the first 
1, 000 iterations have been removed from the trajectory 
to avoid transient behaviour) for different N obtained 
with the same set of thresholds e as in Fig. [7] Dashed 
horizontal lines indicate numerical estimates of Z?^' 
obtained with N — 30, 000 data points. 



For the generalised baker's map, detailed analyti- 
cal expressions are available for the global Lyapunov 
dimension of the system [53 55 1. For the local version 
of this measure, we obtain the following results: 



DL{x,y) 




y < a A a > Xa 
y < a A a < Xa 
y>aAl~a>Xi, 
y>aAl~a<Xi,. 



(43) 



These expressions can be used as a benchmark with 
which we can compare the numerical estimates of our 
new measures D'^'l and D'!p\ 

In contrast to the Lyapunov dimensions, simple ex- 
pressions for the local clustering and transitivity di- 
mensions of the generalised baker's map can (unlike 
for its symmetric version) only be obtained for some 
specific cases. For example, concerning the dependence 
on a for fixed Aq = A?, = 1/4 and e = 7/4" (n > 1), the 
transitivity can be calculated as 



T=l- 



(2a2(l-a)2)(a2 + (l-a)2) 



1 + a(l - a)(2a(l - a)(a(l - a) + 1) - 3) ' 

(44) 

which allows deriving a corresponding expression for 
D-J-- In the derivation of the latter expression, we have 
used the fact that each linked triple lies in a small 
band of width 16/4" that is composed of four substrips 
of width 1/4" with relative weights of a^, a{l ~ a), 
a{l - a), and (1 - a)'^, and gaps of width 1/4", 5/4", 
and 1/4", respectively. For other values of e, the tran- 
sitivity might be even smaller, leading to larger values 
of Dr- Hence, the estimate based on Eq. ( |44| has to be 
considered a lower bound for the actual value of the 
upper transitivity dimension Dj-. 
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Fig. 9. Estimates of the clustering dimensions ob- 
tained for one realisation of the standard baker's map 
{a — 1/2 and \a ~ = 1/4) with initial conditions 
(xo, yo) = (0.2, 0.2). All results have been obtained with 
= 30, 000 data points. Initial transients have been 
avoided by removing the first 1, 000 iterations from the 
sample trajectory. (A) Colour-coded representation of 
the local clustering coefficient obtained with e = 0.1. 
(B) Scale-local estimates of the transitivity dimension 
Dt{^) in dependence on the considered e. Horizontal 
lines indicate the true values of these measures de- 
rived analytically. Qualitatively the same results have 
been obtained for other (lower) choices of N (not shown 
here), indicating reasonable convergence properties of 
the proposed estimator In addition, point estimates of 
the local upper (C,D) and lower (E,F) clustering dimen- 
sions D^'\ (estimated with 100 equally spaced values of 
e in [0.005, 0.1]) are shown in dependence on both vari- 
ables X (C,E) and y (D,F). Solid horizontal lines corre- 
spond to the theoretical values, whereas dashed lines 
indicate the median values obtained from the consid- 
ered sample of state vectors. 



Figure lOA shows the a-dependence for different 
notions of dimension. Analytical results for the "classi- 
cal" measures Do, Di and D2 have been taken from [55|. 
Note that due to its definition II56I . the maximum lo- 
cal Lyapunov dimension is bound from above by 
the dimension of the underlying phase space (to = 2, 
see Eq. (43 0, although numerical calculations would 
yield higher values for a < Aq. We emphasise that 
setting some dimension of a given set equal to that 
of the surrounding space (m) whenever its numerical 
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Fig. 10. Dependence of several global measures of di- 
mensionality on the parameters (A) a (Xa = Xb = 1/4) 
and (B) A^ = Xt (a = 1/2) of the generalised baker's 
map (here, one realisation has been considered for 
each parameter combination, initial conditions and re- 
moval of initial transients as in Fig. [9|. The grey line 
in (A) corresponds to the theoretical lower bound of 
Dj- (see text). Note that in (B), all "classical" measures 
(Dq, Di, D2 and D'£^°-^) have equal values. Numerical 
estimates have been obtained with N = 15, 000 (A) 
and N = 25, 000 (B) data points, respectively, using 50 
equidistant values of e e [0.001,0.015] (in this range, 
and Dij- approximate the analytical results well for 
the a = 0.5 case as shown in Fig. |9j3). The results are 
robust for various choices of N. 



value exceeds m avoids pathological behaviour, which 
has been observed in a similar way for the clustering 
and tra nsitivit y dimensions introduced in this paper 
(see Sec. 3.2.3| l. Concerning the numerically estimated 
transitivity dimensions, we observe that coincides 
rather well with the lower analytical bound resulting 
from (44 1, but shows some positive bias. There are two 
possible reasons for this: (i) either the estimate ( |44| l for 
the lower bound of the transitivity T is still too conser- 
vative, or (ii) the numerical values are too large due to 
some overshooting as indicated in Fig. [9p. Comparing 



Z)^' with the other dimension measures, it becomes 



evident that in general, the transitivity dimension is 
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neither an upper nor a lower bound to any of the con- 
sidered classical concepts. 

The dependence of different dimension measures 
on Xa = Xb is shown for a = 1/2 in Fig. [lOb. Note 



that for this specific choice of a, all "classical" dimen- 
sions take the same values, which only depend on = 
Afc. For sufficiently small Xa = Xf, < 0.3, the upper 
and lower transitivity dimensions seem to take sta- 
tionary values, which is in contrast to the other mea- 
sures of dimensionality. For = — > 0.5, the map 
fills the complete two-dimensional unit box, so that all 
dimensions converge to 2 (note that this limit is not 
approached by the lower transitivity dimension due to 
the finite length of the considered realisation). For too 
small Xa = Xb, the numerical behaviour also suggests 
that longer realisations of the system are necessary to 
obtain reasonable results. In general, we conclude that 
for the parameter range within which our results can 
be considered reliably, the numerically estimated tran- 
sitivity dimensions take similar values as the other 
measures and show a similar behaviour if the param- 
eters of the generalised baker's map are varied (with 
the exception of the maximum local Lyapunov dimen- 
sion), which suggests that the new network-based di- 
mensions are reasonably defined. 

In all cases, note that the numerically estimated 
values of show some residual variations superim- 
posed to their general trend. Besides the finite N, this 
is mainly due to the fact that only one specific realisa- 
tion of the system at every set of parameters is used. 
We expect results to further improve if mean values 
taken from ensembles of independent realisations are 
considered. 



4.3 Rossler system 

So far, we have only discussed examples of discrete 
maps. Among the dynamical systems showing complex 
behaviour, there are however many examples that are 
time-continuous rather than discrete. In the following, 
we will discuss as one paradigmatic example the well- 
studied Rossler system 



X = —y — z 
y = x + ay 
z = b + z{x — c) 



(45) 



with the parameters a — 0.2, 6 = 0.2 and c = 5.7. 
For the latter choice, the Rossler system is known to 
have a chaotic attractor. Moreover, there are count- 
ably many unstable periodic orbits (UPOs) of various 
periods, which do not belong to the attractor, but are 
densely embedded in it and support the invariant mea- 
sure. Hence, these UPOs form a subset of the attrac- 
tor's closure, which has measure 0. As a consequence, 
traditional dimension estimates typically characterise 
the properties of the chaotic part, but are not suited 
for describing the properties of the embedded unsta- 
ble, but dynamically invariant periodic structures. In 



contrast, our results obtained for the logistic map sug- 
gest that local transitivity properties of e-recurrence 
networks can be used for identifying at least the least 
repulsive UPOs. In the following, we will further dis- 
cuss this idea and present some numerical results us- 
ing the concept of continuous clustering and transitiv- 
ity dimensions introduced in this paper. 



4.3.1 General considerations 

Generalising the previous considerations concerning 
the su pe rtrack functions of the logistic map (see Sec. 
|3.2.1| and |4.1.2| l to time-continuous dynamical systems, 
it appears a reasonable assumption that the continu- 
ous e-clustering coefficient C(a:; e) is in general a conse- 
quence of the spatial alignment of neighbouring trajec- 
tories in phase space, which is closely related to the ef- 
fective local dimension of the attractor. In this respect, 
we note that continuous systems may be transformed 
into discrete maps by choosing a proper Poincare sec- 
tion. For example, supertrack-like structures in Poin- 
care sections of the Lorenz system can be identified us- 
ing the vertex properties (in particular, degree and lo- 
cal clustering coefficient) of the associated e-recurrence 
networks [18]. 

Besides our specific considerations for maps, we no- 
te that in general, spatial differences in C{x;e) (and, 
hence, D^^ \x)) can be theoretically understood using 
results for random geometric graphs [61|. Recall that 
since individual recurrence points are assumed to be 
separated by sufficiently large time intervals (i.e., so- 
journ points are excluded) [17 |, the actual spatial loca- 
tion of the associated vertices depends on the specific 
sampling of the data. Therefore, an e-recurrence net- 
work can be interpreted as a random geometric graph 
with a certain effective dimension (in our case charac- 
terised by D-y). We note that the latter considerations 
apply both globally and locally, i.e., they also hold for 
arbitrary subgraphs of an e-recurrence network. Since 
for arbitrary geometric graphs, the subgraph proper- 
ties follow from the spatial distribution of vertices, spa- 
tial heterogeneities in this distribution can result in 
(among others) different local transitivity properties 
and, hence, a non-trivial spatial pattern of the point- 
wise (scale-local) clustering dimension. 

We emphasise that a low local dimension (< m) of 
the attractor implies that trajectories cannot (locally) 
exponentially diverge in all directions of the m-dimen- 
sional phase space, but rather become (locally) almost 
parallel in some lower-dimensional subspace. Among 
other cases, the latter behaviour can be considered typ- 
ical in the vicinity of UPOs, where trajectories become 
dynamically trapped near an invariant lower-dimen- 
sional object for a certain finite time [70 |. Since for 
random geometric graphs, it is known [61 1 (and veri- 
fied by our analytical considerations in Sec. [3]l that the 
expected clustering coefficient decreases roughly ex- 
ponentially with increasing spatial dimension of such 
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networks, the hypothesis that C{x; e) takes local max- 
ima close to UPOs appears justified, which translates 
into a low pointwise scale-local clustering dimension 
Dc{x;e). From this perspective, Dc directly relates to 
traditional concepts like pointwise dimensions (which, 
however, would typically characterise the chaotic at- 
tractor rather than the embedded UPOs) and local Lya- 
punov dimensions (which has, however, only been for- 
mally defined for maps so far). Moreover, we note that 
there is a direct link between Lyapunov dimension and 
Lyapunov exponents, which measure the average di- 
vergence rate of neighbouring trajectories and can be 
used for a local attractor characterisation as well (see 
definition). 



4.3.2 Period-3 UPOs 

As an empirical verification of the above consideration, 
we consider the dependence of the local clustering co- 
efficient Ci{e) on the spatial coordinates of a vertex. 
Specifically, we study the distance of vertices from the 
two period-3 UPOs embedded in the chaotic Rossler at- 
tractor, which are particularly well expressed features 
of the system. 



20 



10 



As it follows from Fig. 11 there is a 
clear indication that close to these UPOs, both vertex 
degree and local clustering coefficient show increased 
values (note that due to the three-dimensionality of 
the system, these maxima are not as well expressed 
as in the case of, e.g., the logistic map, particularly for 
some finite e smearing out the spatial signatures of the 
UPOs). At somewhat larger distances from these in- 
variant objects, we find a clear tendency towards smal- 
ler values of both measures indicated by significant 
negative values of the rank-order correlation coefficients 
PS. For the degree, this is clearly a consequence of the 
trapping feature of UPOs [70 1, while according to our 
theoretical considerations, the corresponding result for 
the local clustering coefficient (and, hence, the associ- 
ated local clustering dimension) is caused by the low 
dimensionality of the UPOs in comparison to the chao- 
tic attractor itself 

For UPOs of higher periods (recall that these are 
densely embedded in the chaotic attractors), we how- 
ever find much weaker signatures in the spatial pat- 
tern of C{x) |17|. This implies that the detection of 
high-periodic UPOs (which are typically more repul- 
sive and, hence, characterised by shorter residence ti- 
mes in their direct vicinity than UPOs of lower pe- 
riod) by means of e-recurrence networks probably re- 
quires longer time series (larger N) and lower recur- 
rence thresholds e. We note that in principle, UPOs 
can also be detected by other types of proximity-based 
complex network approaches to time series analysis, 
for example, cycle networks li32J . 

4.3.3 Bifurcation scenario 

The bifurcation scenario of the Rossler system is very 
rich and shows multiple complex bifurcations between 
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Fig. 11. Point estimates of (A) the local clustering co- 
efficient Ci and (B) the minimum Euclidean distance 
di from the period-3 UPOs of the Rossler system (N = 
10,000, p = 0.01). In addition, the dependence of (C) 
degree density and (D) local clustering coefficient on 
are shown for short distances (N = 50, 000, p = 0.01). 
We emphasise that d has been measured here using 
the Euclidean norm, whereas for the generation of the 
e-recurrence networks, the maximum norm has been 
considered. Note that the respective (rank-order) cor- 
relations ps between vertex properties and distance 
from the UPO are significant and of comparable order 
for both measures. 



periodic and chaotic solutions in dependence on its three 
control parameters. Recently, much interest has been 
spent on the investigation of so-called shrimps II74I75I . 
i.e., specific self-similar periodic windows with a com- 
plex shape that appear in certain two-dimensio nal sub- 
spaces of the full parameter space (see Fig.[l2]l [ 761771781 . 
It has been demonstrated that statistical measures ba- 
sed on recurrence quantification analysis as well as e- 
recurrence networks are well suited for automatically 
discriminating between periodic and chaotic dynam- 
ics and, hence, uncover complex bifurcations between 
both types of behaviour ||37l . Within the complex net- 
work approach, transitivity properties have been found 
to be among the most suitable candidate measures for 
this purpose. Given the framework of our considera- 
tions presented in this work, this effect can be theo- 
retically understood since periodic trajectories corre- 
spond to a lower-dimensional dynamics than chaotic 
ones, which is naturally detected by the transitivity di- 
mension. 



5 Summary 

The recently introduced e-recurrence networks have 
a great potential for detecting qualitative changes in 
the dynamics of complex systems, which may corre- 
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Fig. 12. Variation of the scale-local transitivity dimen- 
sion Dri^) obtained for individual realisations of the 
Rossler system with N ^ 5, 000 and p = 0.02 in a 
two-dimensional cross-section (a = b) of the parameter 
space. Periodic windows are characterised by minima 
of Dr{.£) with values close to 1. White points indicate 
parameter combinations for which the numerical algo- 
rithms did not provide feasible results for the consid- 
ered parameters (e.g., indicated an artificial fixed point 
behaviour due to the improper choice of the sampling 
rate). 



spond to nonstationarities, bifurcations, or different lo- 
cal attractor properties. While previous results have 
been mainly obtained numerically, this paper provides 
a theoretical framework for better understanding the 
links between network and attractor properties. In par- 
ticular, we have studied the local and global transi- 
tivity properties of e-recurrence networks, which are 
closely interrelated with the local and global dimen- 
sionality of the studied attractor. This relationship mo- 
tivated the definition of novel measures of dimension- 
ality, the (local) clustering and (global) transitivity di- 
mensions, which can be directly estimated from this 
type of networks. In this spirit, our corresponding re- 
sults demonstrate that e-recurrence networks provide 
an important link between dynamical systems theory 
on the one hand, and graph theory on the other hand. 

We emphasise that many other established mea- 
sures of dimensionality, such as box-counting and Ren- 
yi dimensions, are based on the scaling of local res- 
idence probabilities of typical trajectories on the at- 
tractor in different parts of the phase space with suc- 
cessively refined spatial resolution. As an exception, 
the correlation dimension is based on spatial two-point 
correlations in phase space. In this respect, cluster- 
ing and transitivity dimensions are statistical proper- 
ties of higher order, since they are based on geometric 
three-point interdependences, i.e., the mutual proxim- 
ity of triples of state vectors on the attractor As a re- 



sult, local clustering dimensions allow quantifying the 
effective (possibly non-integer) local dimensions of the 
attractor in different parts of phase space. The fun- 
damental importance of the corresponding geometric 
interpretation becomes visible in the representation 
of distinct spatial structures related with supertrack 
functions and UPOs, which cannot be detected by other 
traditional measures of dimensionality. 

Beyond the aforementioned conceptual differences, 
we note that our novel dimension measures have fur- 
ther important advantages in comparison to more tra- 
ditional properties such as correlation or pointwise di- 
mensions. These advantages mainly reflect the issue 
of practical estimation: whereas for many classical di- 
mension measures, scaling properties of some quantity 
have to be carefully evaluated (which requires large 
data sets and sophisticated estimation strategies Il62l ). 
there is no need for considering any specific scaling for 
estimating clustering and transitivity dimensions. Be- 
sides the fact that the estimation becomes more direct, 
this also allows numerically obtaining reasonable es- 
timates from rather short time series (i.e., data sets 
of size N ^ 0(10'^ . . . 10"*)) at least for low-dimensional 
systems. The required amount of data is therefore sig- 
nificantly lower than for classical properties such as 
D2, implying that all numerical calculations performed 
for this paper can be completed on standard desktop 
computers within a reasonable amount of time. We ex- 
pect that this advantage of much lower requirements 
with respect to the number of data should persist for 
higher-dimensional systems. 

In contrast to these benefits, we have identified sit- 
uations where our new measures behave pathologically 
(e.g., exceed the non-fractal dimension of the phase 
space in which the attractor is embedded). However, 
we emphasise that similar pathologies may also be found 
for other concepts of fractal dimension, e.g., due to the 
breakdown of the supposed scaling relationships, or in 
terms of the "artificial" upper bound of the Lyapunov 
dimension. The numerical examples discussed in this 
paper demonstrate that there is no simple relationship 
with any existing dimension measure, i.e., clustering 
and transitivity dimensions do not serve as bounds to 
any of the more traditional concepts, but typically have 
values that are comparable with those of other types 
of fractal dimensions estimated from the same trajec- 
tories. 

Our theoretical considerations also confirm recent 
numerical results on the relationship between local tran- 
sitivity properties and the location of dynamically in- 
variant objects. Specifically, for the logistic map, high 
values of the local clustering coefficient coincide with 
the positions of supertrack functions, which has been 
studied in more detail in this work. For the three-di- 
mensional chaotic Rossler oscillator, it has been shown 
that unstable periodic orbits with low periods coincide 
with local maxima of the same vertex property |17). 
Our results suggest that these findings can be gener- 
alised to other (discrete as well as time-continuous) 
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complex systems, given that the invariant density of 
the attractor is sufficiently continuous in phase space. 
Examples such as Cantor sets or the two-dimensional 
Henon map have been discussed as well, illustrating 
the fact that in particular the proper estimation of lo- 
cal (pointwise) dimension measures is non-trivial for 
attractors with a pronounced fractal structure. Our find- 
ings suggest a fundamental relationship between the 
differences of upper and lower clustering/transitivity 
dimensions (which have been found for certain self- 
similar sets) on the one hand, and the smoothness prop- 
erties of the attractor on the other hand. A more de- 
tailed investigation of the corresponding interdepen- 
dences will be subject of future studies. 

The relationship between local transitivity proper- 
ties and local attractor geometry theoretically justi- 
fied in this paper has some important consequences 
for possible practical applications of e-recurrence net- 
works in dynamical systems research. In particular, 
the fact that the local clustering dimensions are excel- 
lent candidates for quantitatively characterising the 
(mean) dimensionality of the system within some e- 
ball around any specific point on the attractor can help 
numerically identifying dynamically invariant objects 
such as unstable periodic orbits (or invariant mani- 
folds of hyperbolic fixed points), which is still a prob- 
lem of intensive scientific research [79|. As a conse- 
quence, we emphasise that our transitivity-based di- 
mension concept offer a novel approach for studying 
structures in the phase space of complex systems and 
appear to have meaningful and potentially relevant 
applications in both dynamical systems theory and real- 
world time series analysis. 
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Table 1. Analj^ical expressions for the expectation values of the local degree densities and clustering coeffi- 
cients for the logistic map at a = 4 in dependence on both x and e [|17i . Note that the remaining integrals cannot 
be solved analytically. 

Starting from the general identities 
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These expressions hold generally for one-dimensional maps defined on the unit interval. 
For the logistic map at a = 4 with the invariant density p{x) according to Eq. ( |30| , one specifically finds 

Ji (0, a; + e) = - — i arcsin(l — 2a; — 2e), 

2 TT 

— £, X + e) = — [arcsin(l — 2a: + 2£) — arcsin(l — 2a; — 2£)] , 



I\{x — e,l) = \ + — arcsin(l — 2a; + 2£), 

2 TT 



72(0,x;0,i/ + £) = 



5- arcsinfl — 2x) , 

4 27r ' ^^J, " ^/{y(l-y)) 



1 arcsin(l - 2j/ - 2£) 



/2(a;, £; 0, a; + e) = ( — \- arcsin(l — 2x — 2£) 1 (arcsin(l — 2a;) — arcsinfl — 2£)) . 

27r TT^ 



/2(£, a; -I- £; y — £, a; + e) = 

72(x — £, a;; a; — E, y + e) = 

72(x, a; + e; y — £, a; + e) = 

/2(a; — £, 1 — e; a; — E, y + e) = 

12(1 — £, a;; a; — £, 1) = 



1 



2a; - 


2£) 


— arcsin(l 


-2e)) + 


^2 




2a; + 


2e) 


— arcsin(l 


-2a;))- 




/ dy 

J X — £ 


2a; - 


2e) 


— arcsin(l 


-2a;)) + 


j^2 


rx + e 

/ dy 

J X 


2a; + 


2e) 


+ arcsin(l 


-2e))- 


j[-2 


r ^ dy 

1 X — E 



arcsin(l — 2y + 2e) 
arcsin(l — 2j/ — 2e) 
arcsin(l — 2y + 2e) 
arcsin(l — 2y — 2e) 



1 ^ 1 

2tv TV 



H arcsin(l — 2x + 2e) J (arcsin(l — 2a;) + arcsin(l — 2e)) . 



hix, 1; 2^ ~ ^' 1) = ^ + ^ arcsin(l — 2a;) + 



dy 



arcsin(l ~ 2y + 2e) 



